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Abstract 

We show both experimentally and numerically that the time scales separation introduced by 
long range activation can induce oscillations and excitability in nonequilibrium reaction-diffusion 
systems that would otherwise only exhibit bistability. Namely, we show that in the Chlorite- 
Tetrathionate reaction, where the autocalytic species diffuses faster than the substrates, the 
spatial bistability domain in the nonequilibrium phase diagram is extended with oscillatory and 
excitability domains. A simple model and a more realistic model qualitatively account for the 
observed dynamical behavior. The latter model provides quantitative agreement with the experi- 
ments. 
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Dynamics of spatial and spatio-temporal react ion- diffusion patterns in chemical systems 
kept far from equilibrium by a permanent feed of fresh reactants is an active domain of 
nonlinear science |]1], 0, ^. These patterns include travelling waves in oscillating or excitable 
media, Turing patterns, self-replicating spots, labyrinthine patterns and dynamically struc- 
tured fronts. They result from the destabilisation of a trivial state of high symmetry by a 
positive feedback ('activation'), caused by autocatalysis or substrate inhibition, competing 
with a negative feedback on appropriate time scales. Although some patterns are transiently 
observed in some systems, these dissipative structures can only be sustained in permanently 
fed open reactors. Among these, the One Side Fed Reactor (OSFR) made of a thin slab of a 
chemically inert gel fed on one face and described further on in detail has become popular. 
Recently, a new type of dynamical behavior, namely spatial bistability, has been evidenced 
in this type of reactor |Q. 

Up to now, two classes of spatial patterns have effectively been observed. In the first class, 
the primary source of the patterns is the dynamical properties of the reaction in homogeneous 
conditions. Reactions which oscillate or are excitable in homogeneous conditions, i.e. in well 
stirred reactors, give rise to travelling waves in unstirred media. The role of diffusion is to 
convert the temporal dynamics into spatio-temporal dynamics 0. The phenomenon may 
occur whatever are the relative values of the diffusion coefficients of the different species. In 
the second class, the source of patterns is the so-called 'long range inhibition'. An activatory 
species, i.e. a species which controls the positive feedback, hke an autocatalytic species, 
diffuses much slower than the inhibitory species. This favors the growth of the instability 
at localized positions by creating two different time scales for spatial interactions The 
prototypes of such patterns are the stationary Turing patterns [0, ||, |^, but long range 
inhibition is also involved in other phenomena like self replicating spots [O, labyrianthine 



patterns |jTT| or morphological instabilities of fronts [|T2|, |T3], |TJ]. Since all small molecules 
in solution have almost the same diffusion coefficients, it was necessary to slow down the 
effective diffusion of the activatory species by complexation with immobile sites to meet 
these conditions experimentally [0, [1^ . 



Actually, 'long range activation' where an activating species diffuses faster than an in- 
hibitory species defines a third class of systems with several diffusive time scales. Like with 
long range inhibition, one expects this differential diffusion be responsible for new insta- 
bilities that are not present in the sole chemical kinetics. However, rather than inducing 
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the growth of the perturbations more or less locally, the fast diffusing activator tends to 
excite perturbations at large distance from the initial inhomogeneities and to transform a 
stationary state into an excitable or an oscillatory one. In a recent study of spatial bista- 
bility in the chlorite-tetrathionate reaction, which exhibits long range activation, we have 
reported a preliminary observation of such a behavior |]T^. Here, we report more extensive 
experimental data and detailed numerical studies of these phenomena. 

In a first part, we summarize former basic results on spatial bistability in OSFRs, and 
present the reaction, the models and the experimental systems. In a second part, we report 
the results obtained with a simplified model to illustrate the general ideas. In a third part, 
we report both experimental results and numerical simulations with a more complete model. 
The last part is devoted to the discussion of the results. 



I. THEORETICAL AND EXPERIMENTAL MATERIALS 



A. Spatial bistability 

An OSFR is made of a film of a chemically resistant gel, of uniform thickness /. One 
face is kept in contact with the homogeneous contents of a continuous stirred tank reactor 
(CSTR). The other face is pressed against an impermeable wall (Fig. |l]). The coupled 
dynamical equations for the concentrations in the CSTR and in the gel are respectively : 

and 

^ = /.(c) + AV^Q (2) 

where Cjo, Cjh and q are the concentrations of species i respectively in the input flow, in the 
CSTR, and inside the gel, is the corresponding diffusion coefficient, r the residence time 
of the reactor , pv the ratio of the volume of the gel to the volume of the CSTR, and r the 
distance to the CSTR/gel interface. The /j's are the reaction rates. In the right hand side 
of Eq. ([1|), the three terms represent the change of the concentrations per time unit. The 
first term gives the contribution of the reaction. The second one represents the input and 
output flows of the species. It contains all the expandable control parameters of the system, 
namely r and the Cjo- The third term results from the diffusive flux of the species at the 
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interface between the gel and the CSTR and represents the feedback of the gel contents on 
the CSTR dynamics. When the volume of the CSTR is large in regard to the volume of the 
gel {py « 1), this last term can generally be dropped so that the chemical state of the 
CSTR is independent of the state of the gel and the concentrations in the CSTR define a 
Dirichlet boundary condition for Eq. at r = 0, whereas a no-flux boundary condition is 
applied at r = /. To avoid ambiguities, let us make more precise the terms we use. Whereas 
in usual OSFRs the thickness / is small in regard to the two other dimensions, in this paper, 
we want to investigate the behavior of the system along the direction orthogonal to the faces 
in contact with the CSTR and the impermeable wall. Thus, we shall consider systems made 
of 'fiat' gels in which one of the dimensions parallel to these faces {h in Fig. P is still smaller 
and negligible. Such systems can be considered as two-dimensional. Nevertheless, the size 
Z, in the direction orthogonal to the faces, remains small in regard to the second direction 
parallel to these faces. Size / and the direction orthogonal to the faces will be further 
referred to as the 'depth' or the 'width' of the gel, the negligible dimension h becoming the 
'thickness'. 

A number of reactions which present autocatalytic steps are bistable in a CSTR 2, 18 



These reactions normally exhibit 'clock' behavior in batch conditions, i.e. a more or less 
long induction period followed by a fast switch to the equilibrium. In these reactions, 
the transformation rate of the substrate(s) tremendously increases during the switch which 
occurs when the level of the autocatalytic species is large enough. For residence times much 
shorter than the typical reaction time, i.e. at high flow rates, the composition inside the 
CSTR is close to the composition in the input flow. It is referred as the 'Flow state' or state 
F. For residence times much larger than the typical reaction time, i.e. at low flow rates, 
the reaction is almost completed in the CSTR and the composition in the CSTR is close 
to a thermodynamic equilibrium. This is referred as the 'Thermodynamic state' or state T. 
For intermediate flow rates, both states are stable for the same composition of the input 
flows. Spatial bistability can be observed when such a reaction is operated in an OSFR, 
the CSTR of which is set in the flow state. Although this phenomenon, which extends the 



concept of bistability to space, has been previously described in simple models |]T8[ and in 
real experiments [0, its importance has been recognized only recently and a detailed 
analysis was worked out in relation with the development of the OSFRs ^ P^ j. 

To understand the basic principles of spatial bistability, let us consider the concentration 
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profile of a substrate of concentration X and diffusion coefficient D consumed in tlie reaction 
and let us imagine that the depth I of the gel can be changed in a continuous way. The 
time taken by fresh reactants fed by diffusion at the CSTR boundary to reach the opposite 
boundary is of the order of = l^/D. We start from / = and increase / adiabatically so 
that the system is always in its asymptotic stationary state. Let be the induction time for 
the CSTR composition (in state F). While / is small enough for to ti to hold, the whole 
volume of the gel remains in state F. We shall still refer this state of the gel as a F state. 
When to reaches a value of the order of ti, for a value / = /max, the fast reaction starts at 
the impermeable wall and a steep front connecting the F state to a T state forms. Such a 
state will be referred as a 'mixed state' or FT state. But, the amount of the autocatalytic 
species is now large at the impermeable wall of the gel (r = /) and can diffuse backwards 
switching on the reaction, so that the front can recede to a value I = 6 < /max- Actually, if 
one assumes that the reaction rate is extremely large inside the front, one can neglect the 
rate of the consumption of the substrate elsewhere and the concentration of the latter can 
be approximated by the solution of a diffusion equation between a fixed value at r = (the 
CSTR) and r = 6 (the front). Then, it can easily be shown that: 

(3) 

Q 

where Q is the total amount of substrate converted by time unit in the gel and AX is the 
variation of concentration between the F and the T states [Q. This value is practically 
independent of /, so that one can now decrease / down to a value /min ~ S before the system 
switches back from the FT to the F state. Thus, for /min < / < /max, both concentrations 
profiles F and FT are stable. This defines the spatial bistability. The extent of the bistability 
domain depends on the concentrations in the CSTR - always close to those in the input 
flow - and the residence time, so that various nonequilibrium phase diagrams can be built 
0, In spite of its apparent crudeness, this theory gives accurate results. Extensive 

studies have been performed experimentally and numerically on the Chlorine Dioxide- Iodide 
(CDI) reaction. Up to now, differences in diffusion coefficients did not seem to play any 
role. Actually, whether the system is ruled by equal diffusion coefficients or by long range 
inhibition, it exhibits a standard spatial bistability. We shall see that the situation is fairly 
more complex when long range activation comes into play. 
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B. The Chlorite- Tetrathionate reaction 



The Chlorite- Tetrathionate (C-T) reaction has essentially been studied in closed systems 
for the determination of stoichiometric and kinetic aspects [p3| , 2^] as well as a source for 



propagating cellular fronts |I2|, O, WM. In excess of chlorite, the reaction kinetics is well 



approximated by the following overall balance equation 

7 CIO2 + 2 S4O2- + 6 H2O — > 7 Cr + 8 SO^- + 12 H+ (4) 
and the reaction rate by 

VR = = MC102][S40r][H+]2 (5) 

Since, during the reaction, large pH changes are observed, the fast dissociation equilibria of 
H2O and HSOj have also to be accounted for. This rate law exhibits quadratic autocatalysis 
in [H'*']. Although the validity of this rate law has been determined for a limited range 
of pH (4.7-5.5) p3|, it has been successfully used in much larger domains. Nevertheless, 
the original value k = 10^ M~^s~^ had to be adjusted to fit with experiments. The value 
= 5 X 10^ M^^s^^ has recently been proposed to fit the experimental data of the bistability 
domain of the reaction in a CSTR [0 and will be used in the following. 

It is well known that, whereas most species have a diffusion coefficient close to 
10"^ cm^s^^, the autodiffusion coefficient of the proton H"*" is almost ten times larger. Nev- 
ertheless, accounting for the electroneutrality condition and the presence of the slowest ions, 
the effective diffusion of H"*" is lower Solving reaction-diffusion equations that include 



all ionic effects is a formidable task which has been only performed on simple models P6 
For practical reasons, it is reasonable to use a constant and uniform effective value for each 
diffusion coefficient. We shall assume that all diffusion coefficients are equal to a unique 
effective value Di, except for and 0H~ which have both a higher autodiffusion coeffi- 
cient and are closely related. Let D2 be a common value for these two species. With these 
assumptions, it was shown in ref fl^ that the experimental data can be correctly interpreted 
only if d = D2/D1 ranges slightly above 3. Thus it is necessary to take into account that 
the autocatalytic species diffuses much faster than the inhibitory species, so that the 
reaction-diffusion system enters the class of long range activation. 
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C. The simple model 



In order to make extensive computations, we shall first define a simple model, the prop- 
erties of which are close enough to the original system to draw significant conclusions. Our 
main purpose is that, due to its formal simplicity, this model can be further used as a proto- 
type for long range activation. The model relies on two strong assumptions: a) The kinetics 
is completely described by Eqs. (^, b) The concentrations at the CSTR boundary are 
fixed to their values in the input flow, i.e. we assume a residence time r = (infinite flow 
rate). In these conditions, the reaction-diffusion system is: 



- -7A;[C10^][S40^][H+]' + DiV'[C102 



dt 



^^-2Mc.o,||s,onim=.i,.v^|s,on (e) 

In this paper, we do not consider short transients that follow a change of parameters and 
only consider the asymptotic regime. Eliminating the reaction term between the two first 
equations, one is left with a diffusion equation with a fixed boundary condition at r = 0. In 
the long time limit, the solution is 2[C102"] — 7[S40g^] = 2[C102 ]o — 7[S40g^]o, where index 
holds for the concentrations fixed at the boundary. Then [S40g^] can be expressed as a 
function of [C102^] and can be eliminated from Eqs. (^. One gets the two- variable reduced 
system: 

dii 



(7) 



dv 12 , 9 ,„o 

^ = -{u - Ouv' + dVl,v 

where we have introduced the new time and space scales t' = 2A;([C10^]o)^t, 

and the following normalized variables and parameters: 



r 
u 



[a02]/[a02]o, V = [H+]/[a02]o, d = D^/D^ and e = l - (7/2)[S40t]o/[C102']o. This 



system is formerly identical to the system used in previous studies in closed systems [24 



According to the definition of m, one has m = 1 at the CSTR boundary. Thus, the control 
parameters are Wq, rf, and ^, where Vq is the fixed value of v at the CSTR boundary. Pa- 
rameter ^ represents the normalized excess of chlorite. Provided it is not too large - this is 
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the case in most experiments - it was found to play a minor role in the model and does not 
bring qualitative changes in the results. Thus, it was fixed at = which would correspond 
to exact stoichiometric proportions of chorite and tetrathionate. 

Note that, if d = 1, the third equation could be eliminated in the same way and the 
system would be reduced to a unique variable u. A similar reduction to the same one- 
variable equation could also be done for a homogeneous CSTR. A one variable system can 
be bistable but never oscillating or excitable, since two dynamical time scales at least are 
necessary. This suggests that the C-T reaction can be neither oscillating or excitable in a 
CSTR. More sophisticated models support this assertion, which is confirmed by experiments. 

This model suffers from serious limitations. First, at the CSTR boundary, we use fixed 
concentrations equal to the concentrations in the input fiow. Since, in the fiow state, the 
concentrations should be close to those in the input fiows, this seems to be a reasonable 
approximation. Nevertheless, at finite residence times, the values of Vq could correspond to 
non existent or unstable states of the CSTR. The approximation is strictly valid only in the 
practically unreachable limit of infinite fiow rates (r = 0). In the same way, one can fix vq to 
arbitrary high values for which the CSTR would necessarily be in a T state. This can lead to 
ambiguities in the definition of the spatial states F and T and to experimentally unreachable 
regions in phase diagrams so that some care is necessary for their interpretations. Finally, 
values of Vq can also be taken arbitrary low, down to zero, an unrealistic situation since the 
dissociation equilibrium of water imposes that remains positive and finite. Thus, this 
model must be only thought as a 'toy model', with the advantage that the kinetics remains 
closely related to that of a real system. This will allow us to evidence dynamical phenomena 
induced by long range activation. Such models have also been recently used by Benyaich 
et al. in studies of spatial bistability To make direct comparison with experiments, a 
more elaborated model has to be devised. 

D. The extended model 

To avoid oddities due to the simple model oversimplifications, one has to account, not 
only for the dynamics of the CSTR at r = included in the Eq. but also for the way 
in which the H+ ions are introduced in the reactor and the equilibria in which they are 
involved. In the extended model that we used to make comparisons with the experimental 
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results, the following fast equilibria are taken into account in addition to Eq. (1): 

H+ + OH- ^ H2O 



(8) 



with rate law 

Ve = k:-K [H+][OH-] (9) 

and 

H+ + SO^ ^ HSO4 (10) 

with rate law 

Va = K [HSO4 ] - k- [H+] [SOM (11) 

Discarding the species that do not contribute to kinetics in Eqs. ^ |10D, one retains the six 
variables [H+], [OH"], [C102^], [S40g-], [SO4"] and [HSOj]. The following values are used 
for the kinetic constants: k' = 1.4 x 10^^ M s"\ k+ = k' with = 10"^^, k' = 10^^ 
M s-\ fcj" = Ka k~ with pKa = —\og{Ka) = 1.94. The diffusion coefficient Di is fixed 
to a standard value Di = 10^^ cm^s^^ and the fast species diffusion coefficient D2 is fixed 
to D2 = 3.4 X 10^^ cm^s^^ to fit at best the experimental data. Accurate one-dimensional 
numerical solutions of the resulting reaction-diffusion system have been obtained by an 
implicit finite difference method for stiff systems with adaptative grid in time and space. 
Unfortunately, due to the higher complexity and the extreme stiffness of the model, it was 
impossible to make systematic studies of this model in higher dimensions. Two-dimensional 
computations were limited to a few points in parameter space. In order to avoid problems 
in the definition of the input concentrations and to make easier the comparisons with the 
experimental data and a previous work, we use, as the main control parameter, the quantity 
a defined in the experimental section below. 

E. Experimental conditions 

In order to resolve the concentration profiles in the depth / of the reactor and to minimize 
both technical difficulties and parasitic edge effects, we have developed one-side-fed flat 
annular reactors, i.e. annular OSFRs. Complete descriptions and technical details on such 
reactors can be found in ref. Q]. These OSFRs are made of a flat annular piece of 2% 
agarose gel, of thickness h <0.5 mm. They have a fixed outer radius i?out=12.5 mm and 
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an inner radius -Rin such that the width / = i?out — Rm ranges from 0.5 to 3 mm. The 
agarose gel can be considered as chemically inert in regard to the reaction. In this soft 
hydrogel, the diffusion coefficients are close to those in pure water. The ratio Rin/ 1 is large 
enough for the curvature effects to be negligible. The inner rim of the annulus is tighten 
against the impermeable boundary while the outer rim is in contact with the contents of 
a CSTR of 25 cm'^, fed through a single inlet port with appropriately premixed solutions 
of reagents stored in separated reservoirs. The residence time of the reactor r=600 s, the 
temperature T=25°C, and the input flow concentrations of sodium chlorite [NaC102]o = 

I. 9 X 10~^ M and potassium tetrathionate [K2S406]o=0.5x 10^^ M are kept fixed during all 
the experiments. Besides the two above mentioned reagents, variable mixtures of perchloric 
acid and sodium hydroxide are added to the input flow of reagents to control the pH of 
the premixed feed solution. This is done by pumping different volume ratios of acid and 
base, keeping the overall flow rate constant. For experimental and graphic convenience, this 
relative distribution is characterized by a control parameter a on which the acid and base 
flow concentrations depend linearly: one has [HC104]o = a x 0.67 10^^ M and [NaOHjo 
= (1 — a) X 3.33 10~^ M. The pH of the feed mixture decreases from a basic to an acid state 
as a increases from to 1. The parameter a effectively controls the pH of the solution and 
is the only expandable control parameter in this series of experiments. The patterns are 
made visible by bromophenol blue (1% solution), a pH sensitive indicator that turns from 
dark purple in the basic state to pale yellow in the acid state. Patterns are monitored by a 
CCD camera connected to a time lapse video recorder. 

The dynamics of the C-T reaction operated in the CSTR exhibits bistability in the range 
0.18 < a < 0.84 between a F branch (with pH values typically above 9) and a T branch 
(with pH values typically below 3). The upper limit corresponds to the value of a where 
the input flow switches from a basic to an acid solution. The lower limit was used to adjust 
the rate constant k to the value A; = 5 x 10® M^^s"^ [|l^. 

II. THE SIMPLE MODEL ANALYSIS 

In the preliminary experiments of ref. [|l^, we had observed, in addition to spatial bista- 
bility, oscillations of the FT state front and excitability phenomena in the domain where the 
sole F state is stable. While the stability of concentration proflles can be determined with 
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1-D simulations, the latter needs 2-D computations. Thus, the stability limits are deter- 
mined in 1-D calculations by the same accurate implicit method as the full model and the 
existence of excitability domain is studied in two dimensional simulations, performed with 
two different numerical methods, namely an implicit finite difference hopscotch method with 
automatic step control and a finite element method. 

First, choosing a typical depth value / = 10 (expressed in reduced units precedently 
defined), we have built a stability diagram in the plane {vo,d > 1) (Fig. |). At low 
values of Vq and d, there is a large region of spatial bistability that extends down to Vq = 0. 
When vo — > the 'acid' state must obviously be obtained by a continuous decrease of vq 
from another 'acid' state or by an acid perturbation of the system. In Fig. ^ we have 
represented the two stable concentration profiles through the depth of the gel for a given 
set of control parameters. Contrary to the standard case, in the FT state, these profiles 
are always monotonous and do not exhibit a steep front. The states are nevertheless well 
characterized by the low ('basic') or high ('acid') value of f at r = /, the impermeable 
boundary of the gel. When d ^ 1, the system exhibits temporal oscillations in place of 
bistability. Since the system cannot oscillate in homogeneous conditions, these oscillations 
clearly result from the different time scales introduced by the differential diffusion. Although 
the problem is of spatial nature, the stability limits present the characteristic topology of 
the so-called 'cross-shaped diagram' that connects bistability and oscillations of numerous 
reactions in homogeneous conditions PB| . 

A second diagram in the plane {vq, I) has been extensively computed for d=5.5, a value 
chosen for the richness of the dynamical situations exhibited at this diffusion ratio (Fig. 
For small depths, the system is always in the flow (or 'basic') state, as expected. At inter- 
mediate values of depth, there is a large domain of oscillations. It vanishes at large vq values 
for which one continuously changes from the 'basic' to the 'acid' state when / is increased. 
For large depths, the FT (or 'acid') state is again present down to vq = 0, which defines a 
domain of spatial bistability at small vq values. The most interesting part of this diagram is 
the region labeled 'Bist.2' in Fig. §(b). In this narrow region, the system still exhibits spatial 
bistability, but the 'acid' state becomes oscillating. Thus, the bistability domain is divided 
in two regions, one of standard bistability between two stationary concentration profiles 
('Bist.l'), and a second one between a stationary profile and an oscillating one ('Bist.2'). 
Moreover, when one decreases / and crosses the stability limit of the oscillatory 'acid' state. 
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the unique stable state, namely the F ('basic') state, exhibits particular properties. If, in 
this stable state, an acid perturbation (such as f = 1) is applied close to the impermeable 
wall, the perturbation does not smoothly decay to the stable state but stimulates a further 
increase of ions concentration in its vicinity. A pair of pulses forms and propagates 
along the wall in opposite directions. After relaxation of transients, they travel with con- 
stant shape, amplitude and velocity, which precisely defines the phenomenom of excitability. 
The limits of this excitability domain are reported in Fig. §(a). This domain extends all 
the way to the oscillation zone. In the vicinity of this latter zone, the recovery tail of the 
pulse can become oscillatory in space. Although, for more efficiency, these limits have been 
accurately computed in a rectangular reactor, we have also also performed some runs in 
an annulus, analog to that of the experimental device. An example of an excitation pulse 
propagating in such an annulus (partial view) is represented in Fig. |^. The limits of this 
excitability domain are reported in Fig.|^(a). This domain extends to the oscillation zone. 
In the vicinity of this latter zone, the recovery tail of the pulse can become oscillatory in 
space. 



III. EXTENDED MODEL ANALYSIS AND EXPERIMENTAL RESULTS 



In this section, we report the numerical results that have been obtained with the ex- 
tended model and simultaneously make a direct comparison with the experimental data. 
We consider the existence and stability of the F state and the FT state in the plane (a, /). 

Spatial bistability is observed over a significant domain of the plane. For instance, in 
Fig. 1^, we show the two stable concentration profiles (F and FT states) computed for a 
same set of control parameter corresponding to a point located inside the bistability do- 
main. Contrary to the case of the simple model, the FT state exhibits a steep front. 
The concentration profile of the autocatalytic species is plotted in logarithmic scale to 
emphasize on the position and the stiffness of the front. The concentration of S40g~ and 
OH^ are represented in linear coordinates to underline the linear drop of the substrates 
concentrations as predicted by Eq. (|]). 

In experiments, the pH color indicator reveals the two different states [17|. In the FT 
state, the gel is found to remain in a basic state only in a narrow band located along the 
outer rim. The nature of the different dynamical states were determined unambiguously in 
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the whole parameter range that was explored. 

In Fig. 0, we report the computed limits of the bistability domain (full lines) in the plane 
{a, I) and the experimental points. Gel annuli with / < 0.5 mm could not be technically 
achieved, so that this part of the diagram cannot be explored experimentally. Numerically, 
the upper limit Omax of the FT state is slightly lower than the limit ot = 0.833 computed 
for the F state in the CSTR. Actually, in the FT state, the concentration gradients of the 
substrates at r = are large as seen in Fig. |^, so that, although <^ 1, the last term in 
Eq. (|lD is not completely negligible and the contents of the CSTR are entrained into the T 
state slightly before the limit a = q;t is reached. Experimentally, this minor deviation is too 
small to be detected. There is a first major difference with the standard GDI reaction [Q: 
The spatial bistability persists at large values of /, which seems in contradiction with our 
theoretical presentation asserting that the F state cannot exist when the width of the gel is 
large. Actually, this results from the extreme slowness of the reaction when the system is 
in the basic F state. This leads to a quasi infinite induction time. Thus, one cannot switch 
spontaneously from the F state to the FT state by a continuous change of /, if one remains 
within a realistic range. Nevertheless, one reaches this state by an appropriate perturbation 
or choice of the initial condition in the gel, but the bistability domain extends to huge values 
of the depth /. 

Contrary to the standard case of the GDI reaction, the limit of the FT state when a 
decreases at constant I is complex and exhibits dynamical phenomena both numerically 
and experimentally, provided that / is large enough as seen in Fig. |^. When the value of 
a decreases and crosses the dotted line, the position of the front begins to oscillate and S 
becomes a function of time. This domain of oscillatory dynamics is narrow and difficult 
to explore precisely in experiments. However, one dimensional numerical simulations allow 
for a precise exploration of this oscillatory domain. In Fig. ^ we represent the numerical 
oscillations for decreasing values of a at 1=1 mm, but similar sequences are observed for 
other widths of the gel. When a decreases, the amplitude of the oscillations increases 
and they take a more and more relaxation character. When the a values are drawn still 
closer to the stability of the FT state, one observes, over a very narrow range of a values, a 
period doubling sequence leading to chaotic oscillations. When a decreases from the onset 
of oscillations, the period of the oscillations slowly increases. On the contrary, the period 
strongly increases when the width / increases. Since the experimental conditions are never 
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strictly uniform all along the outer rim of the gel and since the width / cannot be perfectly 
constant, the oscillations of the front lead to the formation of pacemakers that emit travelling 
waves with periods comparable to those found in the computations. 

Beyond the stability limit of the FT front, dynamical phenomena persist in the form of 
excitability. In a large domain of parameters, an acid perturbation, made at a localized place 
of the gel annulus, propagates under the form of two opposite travelling pulses along the 
inner rim. After a short transient, each pulse takes a stationary shape and travels without 
deformation along the annulus (Fig. |]). The length of the pulse tail decreases as a decreases 
until the excitability is eventually lost. When the pulses meet, they anihilate each other, but 
by adding a basic perturbation at the right place and time , it is possible to kill one of the 
two pulses. Then, the surviving pulse can rotate indefinitely along the annulus as previously 



done with the Belousov-Zhabotinski reaction, the paradigm of excitable reactions 30 



This behavior is identical to that observed with the simple model and is also reproduced 
with the extended model for the two different sets of parameters that were checked for 
excitability. 



IV. DISCUSSION AND CONCLUSION 



It is known that, in homogeneous bistable systems, the addition of a negative feedback 
that introduces a new time scale can induce oscillations. In the same way, excitability can 
arise from the coupling of a bistable subsystem that exhibits a S-shaped nulcline with the 
feedback dynamics of a recovery variable. 

Our simple model has clearly shown that the additional time scales introduced by long 
range activation can play a similar role and are able to induce excitability and spatio- 
temporal oscillations in a system which can only be bistable in homogeneous conditions. 
When the diffusion of the activator is not significantly higher than the diffusion of the other 
species, the system can only exhibit spatial bistability. But, when the ratio d becomes 
significantly larger than 1, the system exhibits a domain of oscillations connected to the 
bistable region; the diagram has a typical cross-shaped topology In a narrow band 

within the bistability region, the 'acid' state can also become oscillatory. This is analog to 
the experimental observations and to the results produced by the extended model. 

Moreover, an excitability region develops along the stability limit of the FT state. Actu- 



14 



ally, there is no discontinuity between a stable FT state and an excitable pulse with an infinite 
recovery tail within the F state. One must be aware that the recovery is not controlled by a 
slow kinetic process imbedded in the reaction mechanism as for the Belousov-Zhabotinskii 
reaction, but by the diffusion process that rules the feed from the CSTR boundary. Thus, 
the excitability properties do not only depend on intrinsic quantities, i.e. the reaction rates 
and diffusion coefficients, but also on a geometrical characteristic of the reactor, namely the 
depth of the gel. 

Nevertheless, although the simple model predicts that differential diffusion can induce 
oscillations of the FT state in the spatial bistability region and excitability of the F state, 
the global agreement with experimental data, as well as with the extended model, is not very 
good. Except for some analogies in the left upper part of Fig. |^, the shape of the diagrams 
is different. A first reason is that the regions located on the right side of the bistable and 
the 'basic' domains, including the oscillation pocket, are not physically reachable. The 
corresponding values of vo could not be maintained in the CSTR with realistic residence 
times. The CSTR would switch into a thermodynamic state. By the way, with this model, 
it is difficult to speak of a true FT state since the concentration profiles do not exhibit any 
front. This is why we prefer to use the terms 'basic' or 'acid' rather than F and FT when 
speaking of the simple model. The second major reason is that 0H~ and the autocatalytic 
species H"*" are related by a strong constraint independent of the C-T reaction itself, namely 
the fast equilibrium of water dissociation. This leads to buffering processes and excludes 
that vq takes values close to zero, as permitted by the simple model. Thus this toy model, 
which allows for extensive numerical calculations, is an excellent tool to demonstrate the 
mechanisms of long range activation phenomena, but it should not be used for a precise 
analysis of the properties of the real system. In particular, it would be illusive to try to 
connect quantitatively the model parameters with realistic quantities. 

In the extended model, where the CSTR dynamics is included, a direct comparison 
with experiments is possible. Spatial bistability limits are in good agreement. There is 
nevertheless a point that deserves special attention. Species S40g~ and ClOg are the natural 
substrates, but since the system is fed with a basic solution, OH^ must be considered as 
another substrate in spite of its special relation with H^. Since the diffusion coefficients 
of these substrates are not equal, there is an ambiguity in the choice of the species to be 
considered in the basic theory that leads to Eq. (0). Fig. |^ clearly shows that the relevant 
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species is 0H~. Its concentration drops linearly down to very low values at the exact position 
of the front, whereas the other substrates do not exhibit a linear drop all the way down and 
their concentration is far from minimal at the position of the front. A natural explanation 
is that, although 0H~ diffuses fast, reacts much faster with OH^ than with the other 
substrates. This defines the smallest 6 value and controls the overall process. However, the 
tail of the concentration profiles of these other subtrates plays a role in the anticipation of 
the FT — >F transition that is actually found in the computations. 

Like in the simple model, oscillations induced by the long range activation are observed. 
If one decreases a within the bistability domain, the FT state becomes oscillatory when its 
stability limit is approached. The diagram is in good quantitative agreement with experi- 
mental observations. Since oscillations are restricted to a narrow region of control parame- 
ters, the unavoidable small inhomogeneities along the annulus, make impossible to perform 
a fine analysis of this front dynamics. In the model, the complex oscillations and the period 
doubling sequence observed in the 1-D simulations can be understood using the following 
arguments. Actually, coupled oscillators may present the same complex dynamical behavior. 
Close to the onset of the oscillatory regime, the motion of the front exhibits simple period- 
icity. The amplitude of the oscillations 5max — ^min is small and the front remains located 
close to the CSTR boundary. But when a decreases, approaching the limits of existence of 
the FT state, this amplitude becomes larger and larger. When 5max is of the order of /, the 
impermeable wall acts as a reflector for the diffusive concentration fields. Then, the front 
dynamics couples with itself, with delay and partial damping. This dynamics of coupled 
oscillators explains the observed dynamical complexity . 

In spatial bistability, it is possible to prepare two regions of the reactors in the two 
different states F and FT side by side and to study the dynamics of the F/FT interface. 
In the standard case, the direction of propagation of this interface depends on the control 
parameters and changes within the bistability domain according to the relative stability 
of the two states F and FT as was experimentally shown with the GDI reaction [^, pT| . 
In the C-T reaction, this situation is radically different. Because the fast diffusion of the 
autocatalytic species excites the medium ahead of the interface, the FT state propagates 
always into the F state, i.e. within the whole domain of spatial bistability. Moreover, 
this behavior can persist beyond the stability limit of state FT. The F/FT interface still 
propagates into state F, but, since the FT state is now unstable, the system ultimately 
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returns to the F state. This gives rise to a recovery tail that follows the interface. If the 
ratio of the diffusion of the activator to the diffusion of the main substrates is significantly 
larger than 1 (here d > 3), the propagation of the FT state balances the return to the F 
state at the level of the interface. Thus, the process goes on and a pulse can propagate 
without damping in time. If d is smaller, or if the system is too far from the stability limit 
of the FT state, the return to the stable state dominates and the pulse dies. This explains 
the domain of excitability observed in the experiments and accounted for by the models. 

Despite the good agreement of the extended model results with experimental data, this 
model still presents some serious limitations. A major drawback of the extended model is 
the use of effective diffusion coefficients which neglects the ionic interactions. The main 
effect of such interactions is more likely to give different effective diffusion coefficients in the 
flow state and in the thermodynamic state. As a result, this should change the position of 
the front and the exact limits of the oscillatory and excitability domains, but we believe that 
our main conclusions should hold. From an experimental point on view, a severe limitation 
is the control of parameter uniformity along the annular gel/CSTR boundary which makes 
the front oscillations to appear in the form of travelling waves. Finally, it has been recently 
found experimentally that the kinetic mechanisms summarized by Eq. (^) breaks down when 
the concentration of chlorite or tetrathionate are significantly increased P^ . 

Beyond these first experimental observations, one can expect to control the front motion 
by a continuous change of the geometrical control parameter /. Experiments along these 
lines are presently in progress both numerically and experimentally. 
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FIG. 3: Spatial bistability in the simple model: profiles of u and as a function of depth. 
Dotted lines: F state ; full lines: FT state. Parameters: 1=10, d=3, fo=0.005. 
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FIG. 4: Non equilibrium phase diagram of the simple model at d = 5.5 in the plane {vq, I), a) Large 
view b) Detail of a. 
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FIG. 5: Excitability in the simple model: travelling pulse, d = 5.5, I = 17, vq = 0.003. Concentra- 
tion map of v. from minimum (black) to maximum (white). 
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FIG. 6: Computed concentration profiles through the gel depth. Parameters: i=lmm, a=0.72. 
Dotted lines: F state; Full lines: FT state. 
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FIG. 7: Non equilibrium phase diagram in the plane (a, I). Experimental points: I monostable 

with non excitable F state, A monostable with excitable F state, • bistable with oscillatory FT 
state, ▲ bistable with stationary FT state, * the CSTR switches to the thermodynamic state. 
Computed results: the full lines are the limits of the bistability domain, the dotted line is the limit 
of the oscillory FT state within the bistability region, + marks are points for which the F state 
has been shown to be excitable. 
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FIG. 8: Computed oscillations of the front position S (^=lmm): Oscillation 1: a=0.700 ; Oscillation 
2: a=0.692; Period 2: a=0.6895 ; Period 4: q;=0.6894; Period 8: a=0.6892 ; Chaos: a=0.6885. 
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FIG. 9: Excitation pulse (experimental): an acid pulse (clear color) propagates to the right along 
the inner rim. Z=0.3mm, a=0.48. 
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